Sensitivity to perturbations in a quantum chaotic billiard 
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The Loschmidt echo (LE) measures the ability of a system to return to the initial state after a 
forward quantum evolution followed by a backward perturbed one. It has been conjectured that the 
echo of a classically chaotic system decays exponentially, with a decay rate given by the minimum 
between the width F of the local density of states and the Lyapunov exponent. As the perturbation 
strength is increased one obtains a cross-over between both regimes. These predictions are based 
on situations where the Fermi Golden Rule (FGR) is valid. By considering a paradigmatic fully 
chaotic system, the Bunimovich stadium billiard, with a perturbation in a regime for which the 
FGR manifestly does not work, we find a cross over from F to Lyapunov decay. We find that, 
challenging the analytic interpretation, these conjetures are valid even beyond the expected range. 

PACS numbers: 05.45. +b, 03.65.Sq, 03.20,+i 



Hypersensitivity to initial conditions is the key ingredi- 
ent of classical chaos. In quantum mechanics, its absence 
led to the study of other features that could be associated 
with the chaos of the corresponding classical system. Cel- 
ebrated examples are the Gutzwiller trace formula for the 
quantum spectral density, the description of the spectral 
fluctuations by the random matrix theory and the rela- 
tion of spectral correlations to transport |l],|| . 

In an alternative point of view, Peres || suggested that 
quantum dynamics should distinguish regular and irreg- 
ular classical dynamics if the time evolution of an initial 
state for slightly different Hamiltonians are compared. 
That is, the sensitivity of a quantum system should be 
searched not by changing the initial conditions but rather 
by perturbing the Hamiltonian. The natural quantity for 
this investigation is the ability of the system to return to 
the initial state \4>) after being evolved with a Hamilto- 
nian 7io for a period t followed by an identical period of 
unitary evolution with —Hi = —(Tlo + S). This defines 
the quantum Loschmidt echo (LE) 



M(t) 



| exp[iHit/7i] exp[— iH Q t/h] 



(1) 



The perturbation S can represent the uncontrolled de- 
grees of freedom of an environment. As in classical chaos, 
the LE is related to a 'distance' between a perturbed and 
an unperturbed evolution of the same initial state. 

In recent years new hints were available due to the 
advances in Nuclear Magnetic Resonance. The LE was 
measured in a many-body system of interacting spins Q 
in a range where it is known to have spectral signatures 
of chaos. A striking finding was that, when interaction 
with the environment and residual interactions are very 
weak, the decay of M(t) becomes independent of the per- 
turbation strength. In this situation, it depends on the 
dynamical scales of the systems, i.e. on Ho- While the 
complexity of the experimental system did not allow for 
a derivation of the characteristic time for these specific 
system, Jalabert and Pastawski B studied the LE in a 



one-body classically chaotic Hamiltonian with a pertur- 
bation represented by a long range quenched disordered 
potential. They have showed analytically that M(t) may 
decay exponentially with a rate given by the Lyapunov 
exponent of the classical system. As condition, the per- 
turbation must be quantically strong to produce statis- 
tically unpredictable changes in the quantum phase but 
weak enough to leave the underlying classical dynamics 
undisturbed. 

More recently, Jacquod, Silvestrov and Beenakker || 
predicted a cross-over from a perturbation dependent 
regime to the Lyapunov one. However, this prediction 
is based on the strong assumption that the perturbation 
lives in a FGR regime; i.e. the local density of states 
(LDOS) is a Breit-Wigner distribution whose width V 
varies quadratically with the perturbation strength. In 
this situation, M(t) for a wave packet and the survival 
probability of an unperturbed eigenstate have a decay 
rate given by T. Both observables would describe the 
same physics if the correlation between states forming 
the wave packet could be neglected. 

Our aim is to determine whether the perturbation in- 
dependent Lyapunov regime and the cross-over from a T 
decay are possible in a fully chaotic system with a clear 
semiclassical description where the presence of the per- 
turbation is not described by the FGR. This occurs when 
there are strong correlations that could be related to clas- 
sical structures which prevents a description in terms of 
a random matrix theory. This perturbation is then said 
to be non generic JtJ and the LDOS can be very different 
from the Lorentzian analyzed in Ref. ^| . Our positive an- 
swer in such a case opens the question of a semiclassical 
interpretation for the weak perturbation regime. 

We consider the paradigmatic disymmetrized Buni- 
movich stadium billiard ||. It consists of a free parti- 
cle inside a 2-dimensional planar region whose boundary 
C is shown in Fig. [j] The radius r is taken equal to 
unity and the enclosed area is I + 7r/4. This system 
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not only has a great experimental relevance but 
also it is fully chaotic one by oposition to the system 
considered in Ref. 0. Besides, it can rule out the dif- 
fusive effect of disorder suspected to affect the behavior 
of M(t) in a Lorentz gas jnj. The classical dynamics 
is completely defined once the boundary is given. On 
the other hand, to address the quantum mechanics, it is 
necessary to solve the Helmholtz equation, V 2 M = fc 2 ^ 
with appropriate boundary conditions. fc M is the wave 
number and by setting h — 2m = 1, fc 2 results the en- 
ergy. The most commonly used boundary conditions are 
the Dirichlet (hard walls) and the Neumann (acoustics) 
conditions. However, we are interested in the possibility 
of perturbing the quantum system without breaking the 
orthogonal symmetry and leaving the classical motion 
undisturbed [0. This is possible using more generalized 
boundary conditions: 



m + e*7(g)^(?) = o, 



(2) 



where q is a coordinate along the boundary of the bil- 
liard (see Fig. |l|), and n is the unit vector normal to 
the boundary. g(q) is a real function and £ the parame- 
ter controlling the strength of the perturbation. Dirichlet 
boundary conditions are recovered when £ = while Neu- 
mann conditions are satisfied in the limit £ — * oo. The 
eigenf unctions and eigenenergies for the case £ = are 
readily obtained by using the scaling method Q . 

In order to compute the LE in this system, a relation 
between the eigenvalues and eigenfunctions for different 
values of the parameter £ is needed. Based on a recently 
developed Hamiltonian expansion for deformed billiards 
jL3|, it is easy to show that the eigenvalues and eigen- 
functions for different values of the parameter £ can be 
obtained from the Hamiltonian Ho + £(£) which is ex- 
pressed in the basis of eigenstates at £ = (from now on 
we will call </> M to these eigenstates), 



9^i) -5 — ^~ 

c On on 



dq. 



(3) 



The function g(q) measures the strength of the change 
in the boundary condition along the contour. Within a 
perturbation theory it would represent the direction and 
strength of a distortion of the stadium |H| . Here we use 



9(l) = 



(1+a) sin(q — 1) + 1 



< g < 1, 
I < q < 1 + tt/2 



with a = — l/(2+7r/2) that could be assimilated to a dila- 
tion along the horizontal axis and a contraction along the 
perpendicular one. Notice that the integral above could 
be viewed as an inner product among the wave functions 
defined over C. This relation defines an effective 
Hilbert space in a window Afc w Perimeter /Area [[Hi. 
The cut-off function $ Mt/ = exp [-2 (fc 2 - fc 2 ) 2 /(fc Afc) 2 ] 



restricts the effect of the perturbation to states in this en- 
ergy shell of width B ~ koAk. It allows to deal with a 
basis of finite dimension with wave numbers around the 
mean value fco and restricting to a particular region Ak 
of interest. 

Figure |l| shows the dependence of the energy levels on 
the perturbation. They exhibit many avoided crossings 
as £ is varied. While the energy levels show the typical 
behavior of a general system without constants of motion, 
we also recognize that some small avoided crossings are 
situated along parallel tilted lines. These energies cor- 
respond to the well known "bouncing ball" states which 
are highly localized in momentum. The selected pertur- 
bation does not modify substantially those states. 

Since the LE is a classically motivated quantity, a 
Gaussian wave packet (with a mean value of momentum 
ko and velocity Vq) is a proper semiclassical selection for 
an initial condition. By evaluating its evolution in a sys- 
tem without perturbation (£ = 0) and other with per- 
turbation strength £, we compute the LE (Eq. I) as a 
function of time. At this point one must recognize that 
the choice of a semiclassical initial condition is very rel- 
evant in order to observe the 'Lyapunov' regime [pTpl] . 

While a global exponential decay of M(t) can be clearly 
identified in almost any individual initial condition, the 
fluctuations for a system with ko not too large can in- 
troduce error in the estimation of the rate. Hence, we 
have taken an average over 30 initial states. Fig. (^) (a) 
and (b) show typical sets of curves of M(t) for ko = 50 
and fco = 100 respectively. It can be seen that after a 
transient, M(t) decays exponentially, ~ exp[— i/r^]. For 
£ > £ c ~ 4.5/A; the decay rate becomes independent 
of the perturbation and 1 /r^ w A with A the Lyapunov 
exponent of the classical system p6| in accordance with 
the conjecture. On the other hand, for large times M(t) 
saturates to a finite value w 1/N with N the effec- 
tive dimension of the Hilbert space M . 

According to Ref. || the chaos controlled decay ap- 
pears provided that A > 1/t where £ = v t is the length 
over which the perturbation changes the quantum phase 
(mean free path) which, for a plane wave with wave num- 
ber k and velocity v . For a quenched disorder perturba- 
tion is evaluated from the FGR 101 



,L = ?1 Um y | Efc , fe | 2 1 



n/2 



(n/2? 



(4) 



Ref. Q realized that in the opposite regime of A < 1/t, 
the LE of an eigenstates <^ of Ho is just a survival 
probability an must decay exponentially under the action 
of the perturbation, 



I (0 M |exp[i(W o +£)i/fi] l<^> I" ~ exphV^ 



(5) 



given by Eq. |j for %/B < t < h/A (A the mean level 
spacing) |lj] . The appearance of this FGR behavior re- 
quires that a typical matrix element U ~ (|S MJ , (£)|} typ of 
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the perturbation to be U > A. The Fourier transform of 
Eq. |5| is the LDOS which, although being discrete, would 
present a Lorentzian envelope [|l8| of width T = 1/r. In 
Ref. H it is conjectured that this decay can determine 
the LE decay with more general initial states. This is 
the regime controlled by the non-diagonal terms in the 
semiclassical expansion || . Once the non-diagonal terms 
have decayed, one expects the chaos controlled decay of 
the diagonal ones will survive. This gives a cross-over 
criterion for the decay rate of the LE of l/r^ = min[r, A] 
as the perturbative parameter £ changes. 

The LDOS is shown in Fig. [| for three different per- 
turbation strengths. In contrast to the case of Ref. H our 
distribution is not Lorentzian. This is related to the fact 
that the used perturbation (the function g(q)) does not 
connect all different regions of phase space; for instance, 
the bouncing ball states are practically undisturbed by £ 
determining the non-generic nature of the perturbation. 
In particular, we have evaluated the width T, showing 
the spreading of the unperturbed eigenstates when ex- 
pressed in terms of the new ones. The results show a 
linear dependence of T on £ shown in Fig. 3; that is, we 
obtain T ~ 0.36£fc 2 . Moreover, taking into account that 
A ~ 0.86/c, the critical value £ c for the crossover from the 
r regime to the Lyapunov one is expected at £ c = 2.4/fc 
(remember that from Fig. ||it results £ c ~ 4.5/k). Then, 
for our system, the criterium works with a T given by the 
half width of the LDOS. This is shown in Fig. ^ where 
for perturbation strengths £ < 4.5fc _1 , the LE. decays as 
M(t) = exp ht/r ] l/Ttj, = T/2 for A > T/2. 

These results contrast with the FGR dependence of 
I/t^ oc £ 2 observed for weak perturbations. These are 
the Lorentz gas with a perturbed effective mass jllj , the 
kicked top perturbed by a perpendicular delayed kick Q 
and general chaotic system perturbed by a quenched dis- 
order || where random matrix theory describes JnJ the 
r decay. In this context, the linear dependence of 1/r 
on £ may be considered as a further indicative that the 
physics of the LE decay can be very different from that 
described by Eq. || and that the result of Ref. || has more 
general validity than expected. In the non-perturbativc 
regime, before the Lyapunov exponent takes over, the LE 
decays exponentially with a rate given by the perturba- 
tion dependent width of the LDOS. Another important 
feature is that £ c ~ 4.5/ k — > when k — > oo. This con- 
firms that in the classical limit Eq. ([j]) would decay with 
the Lyapunov regime regardless of the magnitude of S, 
recovering the chaotic hypersensitivity to perturbations. 

In summary, by studying one of the most important 
models in quantum chaos, a fully chaotic billiard system, 
we have shown that, for a wide range of parameters, the 
Loschmidt echo decays exponentially with rate given by 
the Lyapunov exponent of the classical system. More- 
over, we have discussed the onset of this Lyapunov regime 
requires that A > T/2. In the opposite situation, the pres- 
ence of an exponential controlled by T even in absence of 



a generic perturbation described by the FGR, demands 
further studies to fully interpret the detailed mechanism 
controlling this regime. We finally remark that the M(t) 
would behave much differently for intrisically quantum 
initial conditions. For an eigenstate of ft one finds a de- 
cay described by a FGR and it does not show a crossover 
into the Lyapunov decay [ pT[ . In the other quantum 
extreme, an initial state generated from the long time 
evolution of a semiclassical wave packet p0[ ], we find a 
perturbation dependent Gaussian decay plfl . These is- 
sues have begun to receive much attention p2| ] due to 
its strong connection with quantum computing stability, 
decoherence in waves, and quantum-classical transition. 
Furthermore, the dephasing time observed in transport 
experiments in mesoscopic devices shows a perturbation 
independent rate 0]. So far, there is no consensus about 
the physical phenomenon causing it. Since the time scale 
T(f, measured by the LE is a decoherence time and our 
methodology can obviously be adapted to treat the trans- 
port problem [0], our results open a rich field for explo- 
ration: the connection of both time scales. 
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FIG. 1. Spectrum of the desymmetrized stadium billiard 
with mixed boundary conditions controled by the parameter 
£ [Eq. (3)]. The wave numbers fc M (£) run between 49.3 and 
50.7. Inset: Schematic figure of the system. In solid line we 
show the boundary of the stadium billiard where the mixed 
boundary conditions are applied [Eq. (2)]. The coordinate q 
on the boundary is also shown. Dashed lines correspond to 
the symmetries axis with Dirichlet boundary conditions. 




FIG. 2. M(t) for the desymmetrized stadium billiard per- 
turbed by a change in the boundary conditions. The calcula- 
tions is shown in two different energy regions, (a) Corresponds 
to the region around fco = 50. The value of £ is, from the 
top curve to bottom: 0.019, 0.038, 0.057, 0.075, 0.094, 0.11, 
0.13, 0.15 and 0.17. (b) Corresponds to the region around 
fco = 100. The value of £ is, from the top curve to bottom: 
0.0066, 0.0131, 0.020, 0.0262, 0.0327, 0.0393, 0.0458, 0.0524, 
0.0589, 0.066, and 0.072. The thick lines corresponds to an 
exponential decay with decay rate — 1/A. 
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FIG. 3. Width F of the local density of states as a function 
of the perturbation strenght £ for fco = 50 (filled circles) and 
fco = 100 (circles). The solid lines are the best linear fit. 
Inset: Local density of states P(r) for diferent perturbations 
in fco = 50 (r is measured in mean level spacing units). 
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FIG. 4. M(t) as a function of the rescaled time £i for 
fco = 50 and £ = 0.019, 0.038, 0.057 and 0.075. The dot- 
ted line gives the decay M(t) = exp(— Tt/2). 
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